Energy dissipation statistics in the random fuse model 
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We study the statistics of the dissipated energy in the two-dimensional random fuse model for fracture under 
different imposed strain conditions. By means of extensive numerical simulations we compare different ways 
to compute the dissipated energy. In the case of a infinitely slow driving rate (quasi-static model) we find that 
the probability distribution of the released energy shows two different scaling regions separated by a sharp 
energy crossover. At low energies, the probability of having an event of energy E decays as ~ E~^^^, which is 
robust and independent of the energy quantifier used (or lattice type). At high energies fluctuations dominate the 
energy distribution leading to a crossover to a different scaling regime, ~ _E~^'^^, whenever the released energy 
is computed over the whole system. On the contrary, strong finite-size effects are observed if we only consider 
the energy dissipated at microfractures. In a different numerical experiment the quasi-static dynamics condition 
is relaxed, so that the system is driven at finite strain load rates, and we find that the energy distribution decays 
as V^E) ^ E^^ for all the energy range. 

PACS numbers: 46.50.+a, 62.20.Mk 



I. INTRODUCTION 

Acoustic emission (AE) in fracture is an example of a 
broader phenomenon known as 'crackling noise' |1]. A sys- 
tem crackles in response to an external stimulus leading to 
energy dissipation in the form of avalanches of events with 
no characteristic size. Exarnples can be found in volcanic 
rocks ild, crurnnling paper yj], superconductors disor- 
dered magnets jstl , or plastic deformation of materials . The 
release of acoustic energy in fracture experiments is related to 
iiTeversible processes, like creation of microcracks and defor- 
mation. 

AE experiments in stressed materials typically give a 
power-law distribution of the dissipated energy, V{E) ^ 
E~^, as the system is slowly driven towards catastrophic 
failure. This can be seen as a small-scale analog of the 
Gutenberg-Richter law for earthquakes. The exponent f3 is 
found to depend on the material and in the literature one can 
find P fa 1.3 in synthetic plaster |7|], /? w 1.25 in paper js], 
(3 « 2.0 in fiberglass and /? « 1.51 in wood [loi [llli, 
P w 1.5 in cellular glass |12], and j3 « 1.8 in paper peel- 
ing experiments Ii 1 3i] . This spread of values indicates that the 
exponent is not universal and it may depend on the material 
characteristics as well as on the loading conditions. 

In simple discrete models the scale-free response is associ- 
ated with a cascade of local failures. As the external load is 
slowly increased weak elements fail to hold the imposed local 
stress and break. Some internal load redistribution mecha- 
nism, whose details depend on the particular model, increases 
the local stress on the remaining elements and may cause fur- 
ther simultaneous local failures. As a result, the material may 
respond in avalanches of failure events whose size distribution 
is often very broad. 
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For obvious reasons AE is a very useful non-destructive 
tool that, if properly understood, may give crucial informa- 
tion about the damage accumulated in loaded materials. In- 
deed, mean-field calculations and simulations in some simpli- 
fied models indicate that the energy statistics follows a power 
law with an exponent that changes as the catastrophic failure 
point is reached ifTil [Tsll . Events recorded near global fail- 
ure show a power law decay with an exponent distinctively 
smaller than those events recorded far away from the break- 
down point. Interestingly, recent experiments on granite under 
different loading conditions have shown a qualitative agree- 
ment with this theoretical observation fl^. Should this be a 
robust and universal property it could be potentially used as a 
tool to easily diagnose the damage in loaded materials. 

Scaling in fracture is apparently very generic and can also 
be observed in other quantities (see for instance Ref. [17] for 
a recent review). The distribution of quiet times between AE 
events also obeys a power law analog of the Omori law for 
earthquakes occurrence |8, 10, [U [ll [H, HI [H]. Scale- 
invariant behavior is also observed in the height-height corre- 
lations of fracture surfaces 1 17]. The existence of scaling laws 
in experiments suggests an interpretation in terms of critical 
behavior as the global failure point is approached. However, 
results are not conclusive and a straightforward link between 
experiments and models is still lacking. 

From a theoretical point of view fracture physics is a chal- 
lenging open problem. Materials are heterogeneous and dis- 
order details may play a crucial role. However, the existence 
of power law behavior makes it appealing to try to describe 
fracture with very simplified models that include the most rel- 
evant symmetries, dimensionality, disorder features, and load- 
ing conditions. 

There has been a great effort devoted to describe frac- 
ture crack surfaces in a simple way by means of stochastic 
field equations for the dynamics of an advancing front, where 
the crack is modeled as an elastic string driven through a 
disordered medium in the presence of non-local interactions 
||20[ I2TI l22i I23I I24I1 . However, continuous crack line mod- 
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els have failed to fully agree with experiments, seemingly be- 
cause certain important aspects of fracture, like nucleation of 
voids in the advancing front, are not addressed in these contin- 
uum approaches. In this respect, lattice models are expected 
to be more appropriate. Lattice models describe the elastic 
soUd as a network of springs with random failure thresholds 
where the displacement field can be either scalar or vector In 
particular, the random fuse model (RFM) |25] has attracted 
a lot of attention as a minimal model for fracture since it is 
capable of reproducing some essential features of the frac- 
ture phenomenon in a very simple setting. The RFM repre- 
sents a scalar electrical analog of the elasticity problem. This 
correspondence is based on the formal equivalence between 
the mathematical form of the generalized Hooke's law for the 
scalar elastic problem, ai = J2j Qij^j^ ™d 'he equations for 
the electrical problem, 7^ = kijVj, where the stress (cr), 
local elastic modulus (5) and strain (e) can be mapped to the 
current density (7), local conductance (fc), and local potential 
drop {v), respectively. This model goes one step further than 
mean-field models, hke the fiber bundle model, by including 
local distribution of stress between the nodes, while preserv- 
ing the simplicity in the introduction of disorder as a distri- 
bution of thresholds in the system. Much attention has been 
devoted to the RFM in the last two decades, including stud- 
ies concerning the influence of disorder ll26l |27[] . presence of 
multiscalingbehavior in voltage and current ll28ll . roughness 
exponents ||2^|2§,[3^[3ll[32l, damage localization |33|, and 
failure avalanche distributions tH llfl HI HI [Mil HI. 



Interestingly, an experimental realization of the RFM has re- 
cently been studied [391 to investigate the influence of disor- 
der on the I-V characteristic and crack length. 

The dynamics of the RFM takes place in the form of 
avalanches of failure events. Therefore, it is tempting to link 
avalanche activity to the AE observed in experiments. Several 
studies have focused on the statistics of avalanches of failure 
events in the fuse model 1 15, 34, 35, 36] and improved simula- 
tion algorithms have recently allowed to study larger systems 
with much better statistics r30l [37i [ssll . However, numerical 
determinations of the energy exponent (3 under different load- 
ing rates or on different lattices are scarce. 

Lattice models typically predict larger energy exponents 
(3 than those obtained in experiments. A number of reasons 
might be responsible for this discrepancy including the brittle 
nature of the models, the existence of dispersion of acoustic 
waves in experiments, or the lack of an effective time scale 
separation between stress relaxation and loading in real sys- 
tems, which is assumed in quasi-static models. 

Here we study the statistics of the dissipated energy in the 
RFM under different imposed strain conditions. We compare 
seemingly equivalent definitions of the dissipated energy. For 
a infinitely slow driving rate (quasi-static model) we find that 
the PDF of the released energy shows two different scaling 
regions. For low energies, the dissipated energy distribution 
is robust and independent of the lattice geometry. We find 
'P{E) ~ and a simple scaling argument is given to 

explain this behavior However, at high energies the distribu- 
tion exhibits a crossover to a regime where V{E) ^ i?^^-^^, 
whenever the energy losses are evaluated over the whole vol- 



ume of the system. We also study a different estimator of 
the dissipated energy that takes into account the energy losses 
occurring only at microfractures. This microscopic estima- 
tor exhibits strong finite-size effects and fails to show critical 
behavior. In a different setting, when the system is driven at 
finite strain load rates, we find that the energy distribution de- 
cays as ViE) ^ E~^ for all the energy range. 

The paper is organized as follows: In Sec.|II]we define the 
model, in Sec.|III]we discuss dissipated energy, ways to mea- 
sure it and we report on our results for quasi-static and non 
quasi-static dynamics in the fuse model. Finally, we conclude 
with a discussion in Sec.lIVI 



II. THE RANDOM FUSE MODEL 

We study a two-dimensional lattice of fuses with unit con- 
ductivity kj = 1 and a disordered threshold current 7^^ that is 
a quenched random variable from a uniform probability dis- 
tribution in the interval (0, 1). An external voltage is imposed 
between two bus bars placed at the top and the bottom of the 
system and periodic boundary conditions are imposed in the 
horizontal direction. The system can be driven by a slowly 
increasing voltage or current, mimicking the experimental sit- 
uation where either strain or stress load can be imposed, re- 
spectively. At each update, the Kirchhoff voltage equations 
are solved to determine the currents flowing in the lattice. All 
fuses with current exceeding the corresponding local thresh- 
old are blown and, once the entire system is below thresh- 
old, the voltage (current) is again increased. An avalanche of 
failure consists of all the fuses blown between two external 
voltage (current) increments. Each fuse behaves linearly un- 
til the local current 7^ reaches the threshold. Then, the fuse 
burns and irreversibly becomes insulator kj = 0, so that burnt 
fuses are no longer able to carry any current. The current is 
redistributed instantaneously after a fuse is burnt so that cur- 
rent relaxation occurs in time scale that is much faster than 
the driving time scale. In our simulations, voltage driving is 
imposed. The described configuration should be compared 
with experimental setups in mode-I fracture where the strain 
is slowly increased. An experimental comparison of the AE 
and critical behavior of fracture precursors for different load 
features (strain vs. stress loading) and geometries has been 
reported in Ref. |10]. Although imposed strain experiments 
indicate an induced plastic deformation in the final stages, no 
differences were encountered in AE distribution between the 
two different loading conditions. 

The need to solve a large system of linear equations for 
each update implies a high computational cost and limited 
the reachable system size and the statistical sampling in this 
model in the past, when the b est p erformance was achieved by 
conjugate gradient methods [40]. Recently, a new algorithm 
based on rank-one downdate of sparse Cholesky factorizations 
has been introduced i4ll l42ll . which can largely reduce the 
computational cost of the simulations. Here we make use of 
this algorithm to study two-dimensional networks of fuses in 
triangular and diamond lattices. We studied systems of linear 
size ranging from L = 16 to L = 256 and 10'' realizations of 



3 



the disorder. 



III. DISSIPATED ENERGY 

In stressed materials elastic energy is stored due to redistri- 
bution of the external load all across the lattice. Energy dissi- 
pation in the RFM occurs in bursts of breaking events which 
can be be compared with the AE observed in experiments. In 
doing so, one is assuming that the main contribution to AE is 
given by the dissipated elastic energy, while dislocations and 
friction are, in a slowly driven experiment, secondary contri- 
butions to AE. Also, it is worth to keep in mind that in real 
systems one expects that only a fraction of this dissipated en- 
ergy leads to the AE observed, while the remaining losses are 
due to other dissipative mechanisms like dispersion or damp- 
ing of acoustic waves, which are not described by purely elas- 
tic models like the RFM. Several ways to define the dissipated 
energy can be env isag ed as we discuss below. 

Pradhan et. al. 03 found that /3 = 2.7 for the RFM on dia- 
mond lattices under stress loading conditions. They calculated 
the electric power dissipation in the fuse model as the product 
of the voltage drop across the network and the total current 
that flows through it. The power dissipation in the electrical 
model is equivalent to the stored elastic energy in a mechan- 
ical system. This is a global definition that assumes that the 
whole volume of the system contributes to the dissipated en- 
ergy. 

Alternatively, following Salminen et. al. Ill, we can define 
the energy lost during a given avalanche event n as 

En ^ V^AGn ~ V^Sn, (1) 

where AG„ is the change in the elastic modulus due to 
the failure avalanche, s„ is the number of broken bonds 
(avalanche size) of the nth event, and Vn is the correspond- 
ing potential drop between the bus bars (strain imposed in the 
sample). This definition makes use of the global strain imposed 
on the system and can be seen as a coarse-grained or macro- 
scopic measure of the dissipated energy. 

Note that both, Pradhan et. al. |15] and Salminen et. 
al. lUd, definitions of the dissipated energy take into account 
the whole volume of the system. However, there is strong ex- 
perimental evidence indicating that AE is actually a localized 
phenomenon in space and time so that energy release actually 
occurs at microfractures |9, 10], and it is not therefore spread 
across the system. This suggests one should consider other 
ways to define the dissipated energy in the model, in particu- 
lar, it may be interesting to study measures of released energy 
that are directly linked to the bonds involved in avalanches of 
local breaking events. 

In this spirit, we introduce here what we call microscopic 
dissipated energy, which is defined as the sum of the energy 
losses at every element of the system involved in the nth fail- 
ure avalanche. This can be calculated by adding up the energy 
dissipated by each individual broken bond, jj = kjV^, where 
Vj is the local potential drop at bond j. Since fuses break 
right at the threshold we can define the microscopic dissipated 




3 I I \ I \ I I I \ I 

10 100 1000 10000 

t 



FIG. 1: (Color online) Temporal evolution of microscopic energy 
until breakdown for a typical realization of the disorder on diamond 
lattices ranging from L = 16 (top curve) to L = 512 (bottom curve). 
A first region with slope a = 2 is followed by a second region dom- 
inated by fluctuations. 

energy due to the nth failure avalanche as 

En=Y.{lff, (2) 
3 

where the sum runs over each broken bond within the nth 
avalanche. 

We devote the rest of the paper to analyze the dissipated 
energy statistics in the fuse model for the case of infinitesi- 
mal strain (quasi-static model) and also under different finite 
strain rates (non quasi-static model). We shall be comparing 
the numerical results for the different measures of the dissi- 
pated energy discussed above for triangular and diamond lat- 
tices. One would expect that the three definitions give similar 
temporal behavior for the dissipated energy statistics, apart 
from constant factors. Although this is actually the case in the 
low energy range, dissipation statistics differs at high energies 
for different estimators. 



A. Quasi-static case: infinitesimal strain rate 

Let us first focus on infinitesimal driving. For the sake of 
computational convenience, a voltage V — \ \s fixed between 
the two bus bars. The next fuse to burn is determined by 
maXj (7^/7*'') so the external voltage is increased until this 
fuse exactly reaches its threshold. The new current config- 
uration is calculated according to Kirchhoff equations. This 
rearrangement can cause other fuses to overpass their thresh- 
olds without further voltage increase. All the fuses burnt at the 
same external voltage constitute an avalanche. This process is 
repeated until the network becomes disconnected. 
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FIG. 2: (Color online) Temporal evolution of energy dissipation av- 
eraged over 10* realizations in a i = 128 diamond lattice. All 
the three energies, global, macroscopic and microscopic, exhibit the 
same trend. The three curves have been rescaled by its maximum 
for clarity. 



Quasi-static dynamics results in two very well separated 
time scales, i.e. a fast relaxation process and a slow exter- 
nal driving. As occurs in other systems exhibiting a well de- 
fined time scale separation (as for instance in self-organized 
criticality), the natural time scale is given by the slow time 
scale. Therefore, in the following, time refers to the number 
of avalanches occurred. In Fig.[T]we show a typical realization 
of the temporal evolution for the dissipated energy according 
to Eq. (|2|i for different system sizes. One can see that for each 
realization the released energy grows in time as a power-law 
E{t) t^. At later times, fluctuations around the trend in- 
crease (the larger the bigger the system size) as the system 
approaches the total breakdown point. This dynamic behav- 
ior is very robust and independent of the system size or lattice 
type. 

We also found identical temporal behavior for the macro- 
scopic dissipated energy, Eq. ([T]l, and for the dissipated elec- 
tric energy used by Pradham et. al. for both diamond and tri- 
angular lattices. In Fig.|2]we compare the temporal evolution 
of the average dissipated energy on the diamond lattice ac- 
cording to the three definitions. Notably, this behavior is also 
in agreement with that reported in Ref. [43] for a very differ- 
ent dynamical spring model that included acoustic waves. The 
origin of the robust growth law is perhaps more transparent 
in Eq. where the driving potential is increasing linearly 
with time (as we are imposing a quasi-static dynamics). 

The temporal power-law trend gives significant informa- 
tion about the functional form of the dissipated energy statis- 
tics in the REM. The probability density function (PDE) of 
the dissipated energy at time r in a system of lateral size L 
is given by Pt{E,L) = (1/r) dt6[E{t,L) ~ E], where 
S{u) is the Dirac delta distribution. This corresponds to the 



energy distribution to be observed after the first r avalanche 
events. Note that the probability distribution Vt{E, L) is ex- 
pected to be non-stationary, so it should depend explicitly on 
the observation time t. Let us now consider that the dissi- 
pated energy grows in time as a power-law with some expo- 
nent, L) ^ L^°" P + r]{t), where a captures the scaling 
with system size observed in Eig.[T]and r/(i) is a noise term 
representing the random fluctuations around the trend. The 
details of the noise term ri{t) are not known, but one can ar- 
gue they may depend non-trivially on the interplay between 
the evolving currents and the disordered thresholds. Actu- 
ally, as can be readily seen in Eig.[T] fluctuations are strongly 
asymmetric around the average, which immediately implies a 
non-Gaussian distribution of 77, possibly including non-trivial 
correlations. Despite these difficulties one can perform the in- 
tegral in certain limit, up to certain energy cut-off E^ below 
which fluctuations of the energy are negligible. We have 

Vr{E, L) = T-^ I dt SiL-" e + ri{t) -E] = 
Jo 

= r-^L"^" E^^+^^°' X 

X ds5[s°' +r]{L''/°'E^/°'s)/E -1], 0) 

Jo 

and thus keeping only the lowest-order term we aiTive at 

Pr{E, L) - r-ii'"/"i;-i+i/" (4) 

for T / [L'^ E]^^"" ^ 1. This immediately leads to the exis- 
tence of a characteristic energy scale i?x ^ L~°'t" above 
which energy fluctuations dominate the statistics. It is clear 
that the details of the noise statistics (including the distribu- 
tion and temporal coiTelations) would be required to obtain 
the specific mathematical form of the dissipated energy distri- 
bution above the characteristic energy E^ . Eor the REM we 
have an algebraic growth with exponent a « 2 (see Eig.[Tll, so 
we expect to have a energy distribution decaying as ^ E^^l"^ 
for energies E < Ex- 

We are interested here in the distribution statistics after 
complete breakdown is attained. The characteristic time to 
total failure is expected to scale with system size as Tbreak ~ 
L^, where z is the dynamic exponent. Erom Eq. (|4|l the energy 
statistics after failure, V{E, L) = 7'r=Tbrcak(-^i reads 

V{E, L) - (5) 

for energies below a crossover energy E^ ^ L" ^^'^ . 

Eigure [3] shows the probability distribution, ViE., L), with 
statistics collected up to total failure for the above introduced 
global dissipated energy. Two regions can be readily distin- 
guished. The low energy statistics is in excellent agreement 
with a power-law decay ^ E^^^^ over several decades in en- 
ergy. Similar behavior is observed in Eig. |4]for the macro- 
scopic energy measure defined in Eq. ([Til. 

The dependence with system size of the numerical data ob- 
served in Eigs. [3]and|4]can be better characterized by means 
of a finite-size scaling analysis. The behavior of the energy 
distribution suggests the scaling ansatz 

r{E, L) ^ E-P^L-^g{E/Ex), (6) 
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FIG. 3: (Color online) Probability distribution of the global dissi- 
pated energy for different system sizes. The low energy region de- 
cays as a power-law with exponent = 1/2 and shows a crossover 
at Ex ■ The inset shows data collapse according to Eq. The val- 
ues of K and D are in good agreement with the expected relation 
D — OK. Logarithmic binning has been employed. 

where the scaling function Q{u) ^ const for u ^ 1 and 
becomes G{u) ^ u^^-P-^ for u 1. /3a and /3b are the scaling 
exponents of the distribution above and below the crossover, 
respectively. The crossover energy scales with system size as 
E-^ - with some critical exponent D. 

We can now make use of the theoretical relation we derived 
in Eq. (|5]) to prove that the two scaling exponents n and D are 
not independent. Comparing Eqs. (|6]l and (|5]) one obtains that 
the following scaling relations must be fulfilled: 

D = az — a 
K = z — (J /a, 
/3b = l-l/a (7) 

which immediately imply that D = an. Also, according to 
our estimate a — 2 from Fig. [T] we should have /3b = 1/2. 
This reduces the number of free exponents to achieve a good 
data collapse f4^. 

The insets of Figs. ^ and ^ show a data collapse accord- 
ing to Eq. dill with exponents k = 0.75(4), D = 1.55(5) 
and K = 0.75(4), D = 1.50(5), respectively, and the en- 
ergy exponent below the crossover /?b — 1/2. The fit of 
the scaling function for u ^ 1 corresponds to the difference 
/3b — /5a = —2.25(2), and implies that the scaling exponent of 
the energy distribution above the crossover is /3a — 2.75(2), 
identical within error bars for both energy measures. This ex- 
ponent is to be compared with the one calculated by Pradham 
et. al. for the electric power dissipation in the high energy 
region for the diamond lattice in Ref. 1 15], where they report 
/? = 2.7 over two decades of energy. The macroscopic energy 
defined in Eq. ([T]i not only exhibits the same behavior but the 
same exponents in the two regions indicating both definitions 



FIG. 4: (Color online) Probability distribution of macroscopic dissi- 
pated energy, Eq. Q, for different system sizes and the correspond- 
ing data collapse (inset). The same behavior and exponents as in the 
global definition are found 



are completely equivalent. However, as we show below this is 
not the case for the microscopic energy statistics. 

Figure |5] shows the behavior of the microscopic energy de- 
fined in Eq. (|2|i. Recall that this measure is intended to col- 
lect only those contributions to the released energy coming 
from sites participating in the failure avalanche. We find 
that in the low-energy region the distribution also decays as 
V{E^ L) ^ E^^l"^ . However, in this case we observe that the 
probability does not seem to depend significatively on system 
size, K w 0. Correspondingly, D = olk, « and the crossover 
energy Ey, does not vary with system size. The lack of system 
size dependence of the microscopic energy may be related to 
the fact that the macroscopic and global estimators are sensi- 
ble to the whole volume of the system, while the microscopic 
energy is not. The inset of Fig.|5]shows a zoom of the high en- 
ergy region, where strong finite-size effects are demonstrated 
by the variation of the exponent /3a with system size. Data in 
Fig.|5]obviously fail to exhibit finite-size scaling. 

It is worth to stress here that the PDF for all the three energy 
definitions exhibits identical scaling behavior at low energies, 
^ and is robust to changes in system size and lat- 

tice type. This universality arises from the growth law of 
the dissipated energy which is a feature shared by all the three 
definitions for any system size and lattice geometry. However, 
the lack of scaling behavior with system size of the micro- 
scopic energy has a direct impact on the high energy regime 
of its probability distribution. 



B. Non-stationarity and signatures of imminent failure 

The temporal series of the dissipated energy in the fuse 
model are highly non stationary, as can be easily noticed in 
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FIG. 5: (Color online) Probability distribution of microscopic dis- 
sipated energy, Eq. l|2j, for different system sizes. The inset shows 
a zoom of the high-energy region which decays with an exponent 
around P — 2.5 over two decades for the largest system. Data show 
strong finite-size effects and do not obey good scaling behavior. 

Fig. [T] and so the probability density Vt{E, L) depends ex- 
plicitly on the observation time r. This has been claimed to 
be useful to signal the onset of catastrophic failure Il4 lisll . 
with evident practical applications for diagnosing damage in 
loaded materials. 

In 1 15] the power dissipation avalanche distribution for the 
entire breakdown process was compared with that obtained 
only in a very narrow window around breakdown. In order to 
do this, those authors first computed the average over disorder 
samples of the number of fuses (A'^brcak) blown before catas- 
trophic failure and, for every realization, collected statistics 
from events after almost (A'broak) fuses have blown. A possi- 
ble drawback of that procedure is that, since the time required 
to reach total failure largely varies among different disorder 
samples, one is mixing realizations that are very close to com- 
plete failure with others that are, say, half way into it, which 
obey a different statistics. 

Our procedure to obtain the statistics differs significantly 
from that used in Refs. 1 14, 15] and has the advantage that it is 
not affected by this undesired effect. Moreover, in contrast to 
ifTil [Tsll we want to compare here the distribution of released 
energy until breakdown with that obtained when the system is 
at the very beginning of its evolution and how it changes as 
we approach failure. We proceed as follows. For each disor- 
der reaUzation we let the system evolve up to total breakdown, 
which gives the coiTesponding Tbroak for that particular disor- 
der realization. We then compare the collected statistics with 
that observed for that particular disorder realization up to two 
intermediate times, t = Tbrcak/S and Tbrcak/2, that is, with 
the probability density when only the first one eighth and half 
of the failure avalanches are counted, respectively. So that we 
collect statistics from realizations at the same evolution stage. 



FIG. 6: (Color online) Macroscopic energy distribution for differ- 
ent final observation times. The same exponents collapse the dis- 
tribution independently of the observation time but the slope of the 
high-energy tail is shifted to larger values when one measures further 
away from the breakdown instant. The inset shows the energy distri- 
bution for more intermediate observation times for L — 256. It can 
be clearly observed the slope change at breakdown. 



For each observation time, the probability density 
VriEjL) decays as ^ E^^^^ until the crossover energy 
Ex (t, L). In Fig. |6] we plot the behavior of the distribution 
for the macroscopic energy, Eq. ([T]), for different system sizes 
and different observation times, rescaled according to Eq. 
It can be observed that, while the crossover shifts to larger 
energy values as we approach complete breakdown, the scal- 
ing exponent in the second region is conserved as we increase 
the observation times. However, an abrupt change of expo- 
nent is observed only at the breakdown time Tbroak- This ef- 
fect is perhaps better visualized in the inset of Fig. |6] where 
we plot the distribution data (unsealed) for the largest system 
L — 256 and different observation times. Despite we used 
a different measure of the dissipated energy and a different 
way to collect events these results are in agreement with the 
crossover picture between the two limiting behaviors repotted 
in Refs. fli',^^. This indicates that non-stationary effects of 
the dissipated energy temporal signal may actually be useful 
to characterize damage in stressed materials. 



C. Finite driving rate 

The lack of time scale separation in real experiments has 
been suggested [ 17t] as a possible reason for the discrepancy 
with the typical exponents found in numerical simulations of 
quasi-static models. In order to investigate this point we have 
studied the REM under finite driving rates, so that the model 
evolution is no longer quasi-static. We have analyzed both 
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FIG. 7: (Color online) Temporal evolution of microscopic energy in 
an L = 128 diamond lattice for several finite driving rates. On the 
left panel, it can be observed how, as the driving is increased, the 
energy starts to deviate from growth. On the right panels, one can 
see that, for large enough driving, the energy grows exponentially in 
time. Note the linear-log scale in the latter. 
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FIG. 8: (Color online) Macroscopic energy distribution for L — 
16, 32, 64 and 128 diamond lattices under finite strain increments 
dV/V = 10~^%. A single power-law regime with the predicted ex- 
ponent (3 — 1 is observed over the whole range of energies up to a 
finite-size cutoff. 



stress and strain loading conditions and our results were not 
affected by the loading mode we used. For the sake of brevity 
here we only report on our results for the latter. 

Strain is applied on the system by imposing a small poten- 
tial drop between the bus bars in such a way that all the fuses 
are initially below threshold. The voltage is then increased at 
a fixed rate dV/V letting all the fuses over threshold burn, in- 
stead of the slow driving setup studied above. As before, an 
avalanche is defined as all the fuses burnt between two consec- 
utive voltage increments. In this setup, we can still observe an 
effective time scale separation at the early stages of the evolu- 
tion, while potential increments are small. 

Figure|2]summarizes our numerical results for the time evo- 
lution of the microscopic dissipated energy, Eq. (|2|l, at differ- 
ent strain rates in diamond lattices of linear size L = 128. 
We observe that a power-law trend E t" is satisfied with a 
larger exponent a as the strain rate increases. If the strain rate 
is large enough (about dV/V — 0.01%) the evolution is no 
longer described by a power-law, but it becomes exponential 
in time. The same behavior is observed if either the macro- 
scopic or global energy are used. 

According to our simple calculation in Eq. (IHl, we expect 
that, in the limit of exponential growth, a » 1, the PDF of the 
dissipated energy becomes Vt{E) ~ for E E^- We 
observe that the crossover energy Ex diverges as the strain 
rate is increased. In Fig. [T] one can clearly see that when 
we increase the strain rate the energy fluctuations become 
much smaller. In fact, fluctuations are negligible for the whole 
temporal (energy) range for large enough driving rates, when 
the exponential growth sets in (see Fig. |7ll. This means that 
the fall-off tail of the distribution corresponding to energies 
E ^ Ex is completely washed out in the case of large enough 



strain rates. Therefore we expect the dissipated energy proba- 
bility to be 

r{E) ~ E-^ (8) 

if we let the system evolve up to complete breakdown, r = 
Tbrcak, as Well as for any other intermediate times. 

In Fig. |8]we show the PDF of the macroscopic dissipated 
energy under a strain loading rate dV/V — 0.1% measured 
for all events up to total failure in the diamond lattice for 
different system sizes. Identical results are obtained for the 
global energy estimator (not shown). Our numerical results 
are in excellent agreement with the prediction for finite load- 
ing rates in Eq. (|8]l. The probability distribution for the micro- 
scopic energy also scales as for the whole range of ener- 
gies in the case of finite-driving (not shown), but the crossover 
energy shows no dependence with system size, in agreement 
with our above discussed results D ^ 0. 

The conclusion is that the lack of time scale separation 
leads to a significative change in the distribution of the dis- 
sipated energy in the fuse model, as it should be expected. 
From a physical point of view there are strong differences in 
the system dynamics in the case of infinitesimal driving as 
compared with finite driving. If the system is driven at finite 
rates, relaxation to one of the infinitely many metastable con- 
figurations is not reached before a new perturbation acts on the 
system. This gives rise to a highly nonlinear superposition of 
cascades of released energy instead of individual well-defined 
avalanche events. A finite driving rate generically leads to a 
growth of the dissipated energy at a much faster rate than the 
usual quasi-static dynamics, possibly exponential for any fi- 
nite driving rate in large enough systems. In turn, this fast 
growth takes the crossover energy Ex to exponentially large 



8 



values. The result is that the range of energies in which the 
PDF of the dissipated energy is described by Eq. ([8]) becomes 
very large, actually covering the whole range of energies. 

IV. DISCUSSION 

We have studied the RFM for fracture under strain loading 
conditions. We have focused on the dissipated energy in order 
to test the validity of the model to account for the AE statis- 
tics observed in real experiments in loaded materials. Differ- 
ent ways to define the released energy have been discussed, 
including a microscopic quantity that takes into account just 
the energy losses at each broken bond during an avalanche. 

Our results indicate that, for quasi-static dynamics, the dis- 
sipated energy statistics exhibits two very different regions de- 
pending on the energy scale one is looking at. These two scal- 
ing regions are separated by a typical energy ~ L" 
and obey finite-size critical behavior The low-energy re- 
gion, for E < Ex , is well described by a power-law decay 
V{E) ^ E^^l"^, which is robust and independent of lattice 
geometry. We gave a simple scaling argument showing that 
this robustness is linked to the generic growth law of the 
dissipated energy; a feature shared by all the energy estima- 
tors we studied on any system size and lattice type and that it 
is directly linked to the quasi-static nature of the model. The 
statistics above the typical energy E > E^ crosses over to 
V{E) ~ iJ^^-^^ and ranges several decades in energy. 

Apart from scaling factors, the three energy definitions used 
here were expected to show similar statistics. However, while 
the behavior of macroscopic and global energies can be cap- 
tured by the same scaling exponents and the high-energy re- 
gion exponent is well defined, this is not the case for the 
microscopic energy that, although it obeys the same scaling 
form, shows no system size dependence and a different (size- 
dependent) exponent for the high energy region is obtained. 



Regarding the microscopic energy we introduced here one 
must admit that scaling in the high-energy region is not sat- 
isfying in general terms. Not only the numerical value of the 
exponent depends on the lattice size and fails to exhibit finite- 
size scaling, but also the scaling region covers a very narrow 
energy range. This should be particularly relevant when com- 
paring with real fracture experiments that have shown that dis- 
sipated energy participating in AE is not released all across 
the sam ple, but, quite the opposite, localized at microfrac- 
tures I^JTd]. 

Finally, we also studied the fuse model at finite driving 
rates. It is an often expressed belief that relaxing the quasi- 
static condition might lead to (3 exponents that compare bet- 
ter with experiments. We showed that under finite driving 
the cut-off energy diverges exponentially, so that the scaling 
VlE) E~^ dominates all the energy range at any given 
time for large enough driving rates (possibly for any finite 
driving rate in large enough systems). The conclusion is that 
relaxing the quasi-static condition cannot give account of ex- 
periments, where f3 typically ranges between 1.2 - 2.0 de- 
pending on the material. 

The evidence we have up to now about the RFM indicates 
that it might well be the case that other essential aspects to 
quantitatively account for AE energy exponents in real mate- 
rials, like plasticity, dislocations, friction, damping of acoustic 
waves, etc., are missing in the admittedly oversimplified fuse 
model. 
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